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As a theoretical component of the joint effort with the laboratory of Dr. 
Lou Allamandola to search for potential candidates for interstellar organic 
carbon compound that are responsible for the visible diffuse interstellar 
absorption bands (DIBs), we have performed quantum mechanical 
calculations on the electron absorption spectra of selected polycyclic 
aromatic hydrocarbons (PAH) and derivatives. . This work has resulted in 
one publication in press in Chemical Physics (preprint enclosed). 

PAHs have been postulated as important constituents of interstellar dust 
that connects small carbon clusters and amorphous carbon particles. 
These molecules are thought to be responsible for a variety of spectral 
features observed in the ultraviolet, visible, and infrared regions, including 
DIBs. The validity of this hypothesis, however, requires the knowledge of 
detailed one to one correspondence between the specific PAH species and its 
spectrum. Theoretical calculations carried out on naphthalene and 
derivatives by us provided an important bench mark on a more general 
characterization of such correspondence. 

In the completed project, we have studied 15 different species of 
naphthalene, its hydrogen abstraction and addition derivatives and 
corresponding cations and anions (see enclosed preprint). Using 
semiempirical quantum mechanical method INDO/S, the ground 
electronic state of each species was evaluated with restricted Hartree-Fock 
scheme and limited configuration interaction. The lowest energy spin state 
for each species was used for electron absorption calculations. Results 
indicate that these calculations are accurate enough to reproduce the 
spectra of naphthalene cation and anion observed in neon matrix 
(Samalma F. & Allamandola, L. J. Astrophys. J. 1992, 395, 301). The 
spectral pattern of the hydrogen abstraction and addition derivatives 
predicted based on these results indicate that the electron configuration of 
the n orbitals of these species is the dominant determinant. A combined list 
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of 19 absorptions calculated from 4500A to 10400A were compiled and 
suggested as potential candidates that are relevant for the DIBs 
absorptions. 

Continued studies on pyrene and derivatives have revealed the ground 
state symmetries and multiplicites of its neutral, anionic, and cationic 
species. Spectral calculations show that the cation (^Bgg) and the anion 
( 2 A u ) are more likely to have low energy absorptions in the regions between 
lOkK and 20kK (see attached Figure), similar to naphthalene. These 
absorptions, together with those to be determined from the hydrogen 
abstraction and addition derivatives of pyrene, can be used to provide 
additional candidates and suggest experimental work in the search for 
interstellar compounds that are responsible for DIBs. 
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Abstract In order to preselect possible candidates for the origin of diffuse 
interstellar bands observed, semiempirical quantum mechanical method INDO/S 
was applied to the optical spectra of neutral, cationic, and anionic states of 
naphthalene and its hydrogen abstraction and addition derivatives. Comparison 
with experiment shows that the spectra of naphthalene and its ions were reliably 
predicted. The configuration interaction calculations with single-electron 
excitations provided reasonable excited state wave functions compared to ab initial 
calculations that included higher excitations. The degree of similarity of the 
predicted spectra of the hydrogen abstraction and derivatives to those of 
naphthalene and ions depends largely on the similarity of the it electron 
configurations. For the hydrogen addition derivatives, very little resemblance of 
the predicted spectra to naphthalene was found because the disruption of the 
aromatic conjugation system. The relevance of these calculations to astrophysical 
issues is discussed within the context of these polycyclic aromatic hydrocarbon 
models. Comparing the calculated electronic energies to the diffuse interstellar 
bands (DIBs), a list of possible candidates of naphthalene derivatives is established 
which provides selected candidates for a definitive test through laboratory studies. 
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Introduction 

Polycyclic aromatic hydrocarbons (PAHs) have been postulated as important 
interstellar constituents both in free form, where they are the most abundant 
interstellar organic molecules known, and as the building blocks of interstellar dust 
I- 4 . These molecules in their ionized and/or neutral forms are thought to be the 
origin of a variety of interstellar infrared emission features observed in the 700 cm -1 
to 3100 cm -1 range. The main hypothesis in the model dealing with the infrared 
emission is that PAHs probably present as a mixture of positive ions and some 
radicals as well as neutral species are responsible for most of the discrete infrared 
emission bands through ultraviolet-pumped infrared fluorescence in a collision- 
free environment. 

Moreover, the ubiquitous presence of ionized PAHs in the interstellar 
medium makes them very attractive candidates for the diffuse interstellar 
absorption bands (DIBs) 5 . DIBs number now more than a hundred, and extend 
from about 23000 cm'^ (4400A) to the near infrared 6 . 

While PAHs have been postulated as an important interstellar constituent on 
the basis of the IR emission spectra, the ultraviolet, visible, and near-infrared 
absorption properties of the free, unperturbed species, which are required to 
understand the radiation and energy balance in space, are not known. This has 
motivated our semi-empirical quantum mechanical calculations of the electronic 
structure and spectra of the smallest member of PAHs, i.e. naphthalene and its 
derivatives. 

Numerous spectroscopic studies on naphthalene clusters and its radical ions 
have been carried out in salt media 7 * 9 . However, the electronic spectra of the 
trapped species were strongly perturbed by the ionic solvent. A recent study 
focuses on the spectroscopy of naphthalene and its ions truly isolated in low- 
temperature inert gas matrices where the perturbation of the solvent to the 
species under study is minimum and the environment more relevant to 
astrophysical applications 10 , 11 . Theoretical studies of neutral naphthalene 12 
were facilitated by the well-founded assumption that the entire spectrum is 
derived from excitations within the n electron system. Procedures such as the 
PPP and CNDO/S methods resulted in predictions that were in good agreement 
with experiments. These studies also predicted most of the absorption bands of 
the naphthalene cation 13 - 14 and anion 15 * 20 . However, neither experimental nor 
theoretical study of the hydrogen abstraction and addition derivatives have been 
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reported. 

Thus, we have calculated the electronic structure and spectra of these 
naphthalene radical species in their neutral, cationic and anionic forms. In 
order to validate the methods used, we have repeated calculations for naphthalene 
and its anion and cation, and compared them with experimental results and 
previous theoretical calculations. 

In a recent study 11 , it was found that a small number of DIBs (six in all) 
coincide with bands of the naphthalene cation and that a mixture of different 
PAHs and their derivatives may well explain some of the DIBs. Thus, the main 
goal of this work was to provide theoretical analyses of the calculated spectra for 
comparisons with available experimental measurements and also to provide 
specific candidates among the various possible isomers for laboratory studies and 
comparison with astronomical observations of the DIBs. 

Methods 

To assess the method used for geometry optimization, various semiempirical 
techniques were used for naphthalene by using the program MOPAC.^l AMI, 22 
which gave the optimized geometry of naphthalene closest to experiment, was 
used for the optimization of neutral naphthalene radicals resulting from H atom 
abstraction or addition. There are four such radicals: C^H^-a*, C 10 H7-(3 # , 
C^qHq"£X and C 10 H 9 -3* (a and 3 denote the positions where the hydrogen atom is 
abstracted or added). C s symmetry was imposed during the optimization. Since 
hydrogen addition destroys the aromaticity of naphthalene n system, vibrational 
analyses were performed on the optimized geometries of CjQHg-a* and c 10 h 9 -3* 
to verify the geometries obtained to be true energy minima. Energy and spectra 
calculations were carried out by using the semiempirical INDO/S method. ^3 

Results 

As shown in Table I, AMI amd PM3 gave bond lengths of naphthalene 
closest to the average of the experimentally determined values, as deduced by 
Pauling. 24 Thus, AMI was used for optimizing the geometries of naphthalene 
derivatives. The optimized geometries of the four neutral radicals are shown in 
Figure 1. The geometries of the hydrogen atom addition derivatives are not 
distorted from planarity and are true energy minima since all vibrational 
frequencies were found to be real. 
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The 7i MO energy levels of naphthalene and its cation and anion are shown in 
Figure 2. In the two radicals resulting from hydrogen abstraction, the ^A' state 
was found to be the ground state, in which the unpaired electron occupies the 
nonbonding sp^ orbital (a) localized at the dehydrogenated carbon atom. The ^A" 
states with the a orbital doubly occupied and the 7tg orbital singly occupied were 
calculated to be 52.9 kcal/mol and 56.2 kcal/mol higher in energy for C^o^7' a anc * 
C10H7-P respectively. In the two radicals resulting from hydrogen atom addition 
the unpaired electron in the ^A" ground state occupies the n 5 orbital. Since there 

are no near degenerate orbitals in these radical species, the ^A" states were 
assumed to be the ground states. 

For the cations and anions of the dehydrogenated species of naphthalene, 
there are three low-lying states, 3 A", 1 A", and *A' formed by different orbital 
occupancies and spin couplings. For the cations, the ^A" and ^A" states have two 
electrons occupying the 715 and a orbitals. The ^A’ state is resulted primarily 
from the configuration with both electrons occupying the 715 orbital mixed in with 
a second configuration with the two electrons occupying the a orbital. 

To determine the ground states of these ionic derivatives, SCF and limited 
CISD (configuration interaction with single and double excitations) calculations 
were carried out. In the SCF calculations, the state energies were calculated by 
converging the molecular orbitals of each state to a single electron configuration. 
For the CISD calculations, a set of MOs that are unbiased to either the singlet or 
the triplet state were first obtained by applying the configuration-averaged 
Hartree-Fock procedure^ over the three low-lying states resulting from two 
electrons distributing in the two orbitals (a and Ttg for the cation and a and xg for 
the anion). CISD calculations were then performed for each of the low-lying 
states. Three reference configurations, I I ,.x^ct^>, and I ..n 9 c^>, where x 

is 715 for the cations and X0 for the anions, were included for all three states. 
Using the set of average MOs and the three reference configurations, excitations 
of singles and doubles were allowed within a set of nine orbitals, including eight x 
(712-719) and one nonbonding (a) orbitals. 

The determination of the the ground states of the cations and anions of the 
H-atom addition derivatives C 10 H 9 -a + , C 10 H 9 -P + , C 10 H 9 -a' and C 10 H 9 -P" was more 
straight forward than the ionic species of the H-abstracted naphthalenes. Since 
the unpaired electron occupies the X5 orbital in the 2 A" ground state of the neutral 
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compound, the most stable anion is formed by the addition of an electron to that 
orbital and the most stable cation by extracting an electron from that orbital, 
leading in both cases to J A' states. 

The electron occupancies and the relative energies of the three low-lying 
states of the hydrogen abstracted ions calculated at both the SCF and CISD levels 
are listed in Table II. The relative energies calculated with CISD, which takes 
account of excited configurations that mix into the single electron configurations, 
are similar to the ROHF results. At both levels, ^A" and ^A' were calculated to 
be the ground states of the cations and the anions, respectively. In summary, the 
electron occupancies for the ground states of the 15 naphthalene derivatives are 
shown in Table III. 

Singles only configuration interaction (CIS) using the INDO approximation 
has been shown to accurately reproduce the optical spectra of neutral 
naphthalene^ and other PAHs.^ In this study, the CIS spectra calculations for 
naphthalene derivatives include excitations within the full valence MO set for 
naphthalene and its anion and cation. For the hydrogen abstraction and addition 
derivatives of naphthalene, excitations within 10 it orbitals and about 15 a orbitals 
are included. The optimized geometries of the neutral species were used for the 
spectral calculations of the neutral as well as the corresponding ionic states. 

The calculated spectra of all compounds are shown in Figures 3a-3o. In 
addition, the spectra of naphthalene and its cation and anion calculated at the 
experimental geometry of naphthalene,^ are compared with experiments and 
previous calculations in Tables IV-VI. The dominant configurations of the 
excited states are also listed in the tables in detail. Tables VII to XVIII show the 
calculated spectra of naphthalene derivatives. Only those transitions with 
<40, 000cm* 1 are included in these tables. 

Discussions of Computational Results 

As shown in Table II, both hydrogen abstracted naphthalene cations are 
found to have a triplet ground state, ^A". The open shell singlet states, ^A ", is 
higher in energy at the ROHF level, in consistent with Hund's rule. Inclusion of 
electron correlations in the CISD calculations in the complete k space increases 
the relative energies of ^A" by about 2 kcal/mol. The closed shell singlet states, 
U’, was found to be the highest in energy at both ROHF and CISD levels. 

By contrast, for both anions of the hydrogen abstraction derivatives, the 


6 



closed shell singlet state, ^A', was calculated to be the ground state because of the 
high orbital energy of kq . At both the ROHF and CISD levels, the triplet state was 
found to be higher in energy. The relative energies of the open shell singlet states 
are the highest and increase from the ROHF level to the CISD level, in parallel 
with the results calculated for the cations (Table II). 

Comparing with experiments, 28, 29 the calculated spectra of neutral 
naphthalene and its ions demonstrate the reliability of the INDO/S method in 
predicting the electronic spectra of PAHs (Table IV). Four allowed transitions 
l^Bg u , 2 ^B 2 u , 2^Bg u , and 2 ^B 2 U , were found for naphthalene, with the calculated 
oscillator strengths in agreement with the experimental data. The calculated 
excitation energy for the first absorption (32.3 kK), llBg u , agree with experiments 
almost exactly. The energy for the second state (33.9 kK), 1^B2 U , also agrees with 
experiment very well, only 2 kK lower than the experimental value taken in the 
Neon matrix. I® For the higher energy states, the agreements are not as good and 
the excitation energies were underestimated by as much as 4 kK. However, 
compared to typical ab initio calculations that ignores correlation of a electrons, ^0 
the current results are significantly better in predicting excitation energies, as 
demonstrated by other semiempirical studies. ^ 

Although only single excitations were included in CIS, the dominant 
configurations of each excited state of naphthalene are similar to those found in 
calculations with higher excitations included,^ indicating similar qualitative 
descriptions of these excited states. However, the contributions of these dominant 
configurations are probably over estimated because of the omission of excitations 
of doubles and higher. In this study, two main excitations, la u — »2bgg and 
2 bi u -» 2 b 2 g, were found for the first excited state, l^Bg^ with about half to half 
out-of-phase combination. The first excitation is from the highest occupied orbital 
to the second lowest unoccupied orbital and the second excitation is from the 
second highest occupied orbital to the lowest unoccupied orbital. As shown in the 
orbital energy diagram (Figure 2), these two excitations require about the same 
energy and thus mix almost equally. The small oscillator strength for l^Bg u and 
large oscillator strength for 2^Bg u are a result of this mixing, as illustrated by 
Hoijtink et al.^® By contrast, the two dominant configurations of 1^B2 U and 
2*B2 u> la u -*2b2g an d 2 b 2 u -» 2 bgg, require different energies and these two states 
have oscillator strengths closer to each other. 
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For the symmetry forbidden transitions, the INDO/S energies are lower than 
the ab initio energies but the dominant configurations are mostly the same, 
similar to those states of allowed transitions (Table IV). One state, l^Bgg 

involving excitations to a c* orbital, was also found near 50 kK. 

Similar accuracies were obtained for the calculated spectra of naphthalene 
cation and anion compared to experiments (Table V and VI). Although the 
calculated oscillator strengths for the cation spectrum show discrepancies with 
the experimental values obtained in the Neon matrix, they are in good 
agreement with the measurements in the freon matrix^® and polarized spectra. 
Compared to the existing theoretical results on the naphthalene cation and anion 
using other semiempirical methods, 13-15 these predicted spectra are slightly 
more accurate and more complete. In addition to the transitions corresponding to 
the four excited states of naphthalene, more lower-energy excitations were found 
for the cation and the anion. The extra bands result from the excitations to the 
half filled orbitals, la u for the cation and 2 b 2 g for the anion. For the anion, two 

weak transitions mostly resulting from n—>a* were also found to be near 30kK and 
42.2kK respectively. The oscillator strengths of these transitions were found to be 
small, in agreement with experiment. 

Good agreement of the spectra calculated for naphthalene and its cation and 
anion with experiment allowed us to extend the INDO/S method to the prediction 
of the spectra for the twelve uncharacterized naphthalene derivatives, the 
hydrogen abstraction and addition compounds and their cationic and anionic 
forms. The calculated spectra, together with those of naphthalene and its cation 
and anion, are plotted in Figures 3a-3o. 

Some general similarities and differences can be found among the calculated 
spectra of these hydrogen abstraction and addition derivatives. Similar spectra 
are found for derivatives with the same n electron configurations. For example, 
the cations of both hydrogen abstracted derivatives, C^QH^-a* and CiqH 7 -P + , have 
similar excited state profiles as naphthalene cation. They all include a weak 
transition to the half filled X 5 orbital near 16kK, and a moderately strong 

transition near 37 kK, 714 — » 7 tg. Even the weak transitions between 20-25 kK are 
present in these three states. This similarity is a result of the same electron 
occupancy in the n space of the three cations and the small perturbation of the 
second unpaired electron occupying the nonbonding a orbital of the ^A" states of 
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Ci 0 H 7 -a + and CiqH 7 -P + to the spectra. No transitions associated with this a 
orbital were found because excitations of the k electrons to this orbital are 
symmetry forbidden. 

The spectra of the neutral species, C 10 H 7 -a* and CjqH^P*, are also similar 
to the spectrum of naphthalene. A strong absorption resulted from the mixture of 
n 4~ >n 6 an< ^ k 5~* k 7> was found near 43 kK for both compounds. A weak 
transition, was also found near 38 kK, slightly higher in energy than the 

corresponding peak of naphthalene. 

In contrast to the close correspondence between the spectra of the cations and 
neutral compounds of the hydrogen abstracted derivatives and those of 
naphthalene, for the anions, C^QH 7 -a' and C2 qH 7 -(5', the calculated spectra are 
not similar to that of the naphthalene anion. For example, no corresponding low- 
energy (<25kK) absorptions were found for these compounds. This is because, in 
the ground states -^A', the extra electron occupies the nonbonding a orbital of the 
carbon atom. As a result, the antibonding n orbitals are left unoccupied, which is 
exactly the same electron configuration of the neutral naphthalene. In fact, the 

spectra of CiQH 7 -a‘ and CjqH 7 -P' are similar in some aspects to that of 
naphthalene. 

The calculated spectra of the hydrogen addition derivatives vary according to 
the position of hydrogen addition. Because the disruption of the aromatic 
conjugation in the n space, they do not resemble naphthalene and its ions. 
Among the different ionization states of these derivatives, the anions tend to have 
stronger absorptions near 20 kK and the cation have stronger transitions near 30 
kK. Only weak transitions were found for the neutral compounds in these two 
regions. 

These calculated spectra of naphthalene derivatives can be useful in 
identifying extra absorptions observed in the laboratory. For example, absorptions 
of the naphthalene anion observed in the 20-25 kK region were first attributed to 
hydrogen adduct impurities based on the observation that these transitions are 
non-polarized.^O’^l However, these two transitions were reassigned to the 
transitions of naphthalene anion itself by Shida and Iwata.16 According to the 
current results, the naphthalene anion is not predicted to have observable 
transitions in this region. Instead, the hydrogen addition derivatives, C^qHq-P* 
and CjoHg-a* for example, show weak and moderate absorption in this region. 
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Relevance to the Diffuse Interstellar Bands (DIBs) Issue 

The calculated electronic transitions of naphthalene and its neutral and 
ionic derivatives relevant to the DIBs, are reported in Table XIX. From the 
astrophysical point of view, these calculations indicate that: 

1 . Neutral naphthalene and its dehydrogenated counterparts (C10H7-01 an ^ 
C10H7-P) do not absorb in the visible and, hence, cannot contribute to the 
DIBs. On the contrary, the neutral species CioHg-a and C10H9-P do show 
electronic absorptions in the visible region, and could contribute to the DIBs. 

2. Naphthalene anion and its derivatives obtained either by abstraction or 
addition of a hydrogen atom all show electronic transitions falling in the 
visible or/and near-infrared regions, implying possible contributions to the 
DIBs. 

3. Naphthalene cation and its dehydrogenated counterparts (CioH-7-a + and 
CioH_ 7P + ) do absorb in the visible region and , hence, could contribute to the 
DIBs. CioHg-a + and Cio^g-p*, however, do not absorb in the visible and 
cannot contribute to the DIBs. 

In all, 19 electronic transitions of naphthalene and its derivatives fall in the 
range of the diffuse interstellar bands designating these species as potential 
candidates for the DIBs, as shown in Table XIX. 

One should, however, bear in mind two important facts when making such 
correlations. First, theoretical calculations provide the relative energies of the 
electronic levels of each molecular species and, hence, these values can be shifted 
from the absolute values measured experimentally and which can be correlated to 
the DIBs. For example, in the case of the naphthalene cation (CioH 8 + )> a 
comparison of the calculated energy levels with the corresponding values 
measured when the molecular ion is isolated in a neon matrix 10 , 11 which 
constitutes the most relevant environment to simulate the conditions reigning in 
the interstellar medium (i.e. the least perturbing medium), indicates shifts of 
about 400 cm- 1 (162A at ~ 6500A) and 1400 cm- 1 (318A at~ 4500A) for the two visible 
transitions falling in the range of the DIBs (see Tables V and XIX). This implies 
that one needs to deal carefully with such results which must be used only as a 
tool to preselect plausible candidates for the DIBs and help focus on the species 
which deserve an experimental study. Second, these calculations provide the 
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electronic energy levels of naphthalene and its derivatives but no information on 
the vibronic structure which is part of the real spectrum, was included. The 
vibronic structure is particularly important here to explain the interstellar 
features as indicated by the case of the naphthalene cation (C 10 Hg+) where all the 
6 visible bands which do fit DIBs belong to the same electronic band system (the 
lowest energy system at 6579A (theory )/674lA (Ne matrix)) 10 - 11 . The vibronic 
structure will be obtained from the laboratory measurement of the spectra of the 
species preselected from the theoretical calculations. 

Finally, it can be seen how important is the role of theoretical calculations to 
help clarify a problem as complex as the identification of the DIBs. Faced with a 
large family of molecules (PAHs) including numerous isomers which can each 
show multiple possibilities for hydrogen-atoms addition or abstraction and need to 
be considered in their neutral and ionized forms, one is faced with an enormous 
(if not impossible) task to test the relevance of each candidate to the DIBs. 
Theoretical calculations of the electronic structures of a vast number of these 
species allow one to preselect potential candidates for laboratory study to be 
performed in an environment relevant to astrophysical applications. Then, and 
only then, can a one-to-one correspondence between the measured spectra of PAH 
derivatives and the DIBs be attempted. 

Conclusions 

Comparison of the calculated spectra with experiment shows that the optical 
spectra of naphthalene and its ions can be reliably predicted by the INDO/S 
method. The degree of similarity of the predicted spectra of the hydrogen 
abstraction and derivatives to those of naphthalene and ions depends largely on 
the similarity of the k electron configurations. For the hydrogen addition 
derivatives, very little resemblance of the predicted spectra to naphthalene was 
found because of the disruption of the aromatic conjugation system. 

The study of the astrophysical implications of this work indicates that 19 
electronic transitions of naphthalene and its derivatives fall in the visible and 
near-infrared range suggesting a possible correlation with the diffuse interstellar 
bands (DIBs). The selected species, reported in Table XIX, are tentatively 
designated as potential candidates for some of the DIBs, since the calculated wave 
lengths are in the same region. However, it is stressed that further experimental 
studies in an astrophysically relevant environment are needed to establish a 
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definitive one-to-one correlation between the vibronic bands of a given PAH species 
and some of the DIBs. 
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Figure Captions 


Figure 1. 

Figure 2. 
Figure 3. 


Optimized Geometries of the Hydrogen Abstraction and Addition 
Derivatives of Naphthalene. 

n MO Energy Diagram of Naphthalene and Its Ions. 

Calculated Electronic Spectra of Neutral and Ionic Species of 
Naphthalene and Its Hydrogen Abstraction and Addition 
Derivatives. 
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Table I. Comparison of Optimized Bond Lengths of Naphthalene 
Using Different Methods 



AMI 

MINDO/3 

r c-c ,$ 

1.416 

1.430 

r C-Ca 

1.373 

1.380 

r Ca-Cp 

1.422 

1.449 

r cp-cp' 

1.419 

1.462 


MNDO 

PM3 

Expt # 

1.429 

1.415 

1.401 

1.382 

1.369 

1.372 

1.439 

1.421 

1.422 

1.435 

1.410 

1.412 


# Ref. 24. 

$. The carbon atoms bridging the two rings are labeled as C and C'. 



Table II. Relative Energies of the Low-Lying States of C 10 H 7 
Formed by H Atom Abstraction from Naphthalene 


-H„ -Hp 


Charge State 

Configuration 

ROHF 

CISD# 

ROHF 

CISD # 

+ 1 X A' 

7 I l^ 7 t 2^3^ 7 t 4^ 7 l 5^°^ 

12.7 

12.5 

12.4 

12.4 

lA" 

71 i^2^ 7 t 3^ 7 t 4^ JC 5 1 

5.3 

7.8 

2.6 

4.0 

CO 

> 


0.0 

0.0 

0.0 

0.0 

-1 1 A' 

n i 2 ^ 2 ^ 2 ^ 2 ^ 2 ^^ 2 

0.0 

0.0 

0.0 

0.0 

lA" 

n i 2 ^ 2 ^ 2 ^ 2 ^ 2 *^ 1 a 1 

8.7 

14.8 

6.1 

7.2 

3 A" 

n i 2 Jt2^ 7 C 3^ Jt 4^ 7 t 5^ 7I 6 * 

5.3 

6.1 

4.2 

4.0 


# CISD included 8 MOs, 8 electrons, and 3 reference configurations: 
I ..Tc^a^>, I ..Jt 2 a 0 >, I ,. 7 r^a 2 >. MOs were converged to an average 
multiplicity. 


Table III. Ground State Electron Occupancies of Naphthalene Derivatives 


Naphthalene 

-Ha 

"Hp 

+H a 

+Hp 


Cation 

2 A n (..K 4 2 K 5 h 
3 A "(. . 7:4^5 1 a 1 ) 
3 AX.K 4 2 n 5 l a 1 ) 
1 A'(..7t4 2 it5°) 

1 A , (..7t 4 2 JC 6 °) 


Neutral 

1 Ag(..Jt4 2 «5 2 ) 

2 A'(..7t4 2 ii5 2 <T 1 ) 

2 A'(..7t4 2 Jt5 2 o 1 ) 

2 A"(..Jt4 2 Jt5 1 ) 

2 A''(..n 4 2 n^) 


Anion 

2 B 2g (..7i 4 2 TC 5 2 7t 6 1 ) 

^•A’(..JC4 2 jc5 2 o 2 ) 

lA'(..7C4 2 7C5 2 0 2 ) 

^A'(..^4 2 7t5 2 ) 

^•A'(..^4 2 Jt5 2 ) 
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Table V. INDO/S-CI Calculated Excited States of Naphthalene Cation (Energy in 1000 cm' 
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PPP, Zahradnik et al, ref 14. 

Neon matrix, Salama and Allamandola, ref 10. 

Freon matrix, Shita, ref 16. Oscillator strengths were converted from extinction coefficients with band widths of 1000cm 
Polarized spectrum, Hiratsuka et al, ref 17. 



Table VI. INDO/S-CI Calculated Excited States of Naphthalene Anion (Energy in 1000 cm' 
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PPP, Heinze and Zimmermann, ref 15. 

MTHF, Shida, ref 9. Oscillator strengths were converted from extinction coefficients with band widths of 1000cm" 1 . 

Solution, Brandes and Gerdes, ref 18. Oscillator strengths were converted from extinction coefficients with band widths of 1000cm" 1 
Glass, Hoijtink and Zandstra, ref 30. Oscillator strengths were converted from extinction coefficients with band widths of 1000cm" 1 


Table VII. 


Calculated Spectrum of c 10 H 7~ a+ 


< 3 A" 


state 

E(kK) 


f 

polar 

3 3 A" 

16. 

.1 

0 . 

.094 

X 

4 3 A" 

21. 

.9 

0 . 

.030 

y 

1 3 A' 

24. 

.9 

0. 

.004 

z 

5 3 A" 

25. 

.0 

0 . 

.034 

y 

6 3 A" 

27, 

.0 

0 . 

.014 

y 

7 3 A" 

30. 

.0 

0, 

,009 

X 

8 3 A" 

31. 

.1 

0 . 

,003 

X 

10 3 A" 

34. 

.3 

0. 

.009 

X 

4 3 A' 

35. 

.3 

0. 

,001 

z 

11 3 A" 

37. 

.2 

0. 

.053 

X 

12 3 A" 

37. 

.7 

0 . 

,509 

X 


excitation 


3a" (pi3) 

- 

5a" (pi5) 

5a" (pi5) 

- 

6a“ (pi6) 

3a" (pi3) 

- 

8a" (pi8) 

4a“ (pi4) 

- 

7a" (pi7) 

19a' (Casp2) 


6a" (pi6) 

19a' (Casp2) 

- 

8a" (pi8) 

2a" (pi2) 

- 

5a" (pi5) 

4a" (pi4) 

* 

7a" (pi7) 

4a" (pi4) 

- 

7a" (pi7) 

2a" (pi2) 

- 

5a" (pi5) 

4a" (pi4) 

- 

6a" (pi6) 

5a" (pi5) 

- 

8a" (pi8) 

3a" (pi3) 

- 

6a" (pi6) 

4a" (pi4) 

- 

6a" (pi6) 

5a" (pi5) 

- 

7a" (pi7) 

4a" (pi4) 


6a" (pi6) 

3a" (pi3) 

- 

19a' (Casp2 ) 

4a" (pi4) 

“ 

8a" (pi8) 

3a" (pi3) 

- 

7a" (pi7) 

4a" (pi4 ) 

- 

6a“ (pi6) 

2a" (pi2 ) 

- 

6a" (pi6) 

4a" (pi4 ) 

- 

6a" (pi6) 

5a" (pi 5) 

- 

7a" (pi7) 


11 


coef f . ^ 


0.89 

0.69 

0.11 

0.10 

0.67 

0.14 

0.52 

0.38 

0.69 

0.10 

0.75 

0.45 

0.20 

0.16 

0.56 

0.29 

0.84 

0.40 

0.25 

0.15 

0.10 

0.58 

0.17 


Table VIII. Calculated Spectrum of C 10 H 7 -a* ( 2 A') 


state 

E(kK) 

f 

polar' n excitation 

coef f 

2 2 A" 

27.3 

0.007 

z 

19a 1 (Casp2) 

- 6a" (pi6) 

0.53 





19a' (Casp2) 

- 8a" (pi8) 

0.16 





5a" (pi5) 

- 19a' (Casp2> 

0.11 

6 2 A' 

31.9 

0.008 

X 

5a” (pi5) 

- 7a" (pi7) 

0.40 





4a" (pi4) 

- 6a’ (pi6) 

0.28 





5a" (pi5) 

- 9a" (pi9) 

0.11 

8 2 A' 

34.2 

0.008 

xy 

4a" (pi4) 

- 6a" (pi6) 

0.44 





5a" (pi5) 

- 9a" (pi9) 

0.18 





5a" (pi5) 

- 7a" (pi7) 

0.14 

9 2 A' 

36.9 

0.153 

y 

5a“ (pi5) 

- 6a" (pi6) 

0.89 





4a“ (pi4) 

- 7a" (pi7) 

0.08 






Table IX. Calculated Spectrum of C-^QH-y-a ( 1 A I ) 


state 

E (kK) 

f 

polar 'n excitation 

coef f 

1 ^" 

9.9 

0.004 

z 

19a* (Casp2) 

_ 

6a" (pi6) 

0.67 





19a' (Casp2) 

- 

9a" (pi9) 

0.13 





19a' (Casp2) 

- 

10a" (pilO) 

0.12 

3^" 

23.5 

0.002 

z 

19a' (Casp2) 

- 

8a" (pi 8) 

0.41 





19a' (Casp2) 

- 

6a" (pi6 ) 

0.28 





19a' (Casp2) 

- 

10a" (pilO) 

0.19 





19a' (Casp2) 

- 

9a" (pi9) 

0.11 

4 X A“ 

31.3 

0.009 

z 

19a' (Casp2) 

- 

8a’ (pi8) 

0.42 





19a' (Casp2) 

- 

10a" (pilO) 

0.25 





19a' (Casp2) 

- 

9a" (pi9) 

0.19 





19a' (Casp2) 

- 

7a" (pi7) 

0.10 

2 X A 1 

30.5 

0.093 

xy 

5a" (pi5) 

- 

6a" (pi6) 

0.71 





5a" (pi5) 

- 

7a" (pi7) 

0.14 

3^' 

31.9 

0.042 

xy 

5a" (pi5) 

- 

7a" (pi7) 

0.47 





4a" (pi4) 

- 

6a" (pi6) 

0.28 





5a" (pi5) 

- 

6a" (pi6 ) 

0.19 

4 2 A' 

36.4 

0.005 

xy 

19a' (Casp2) 

- 

20a 1 

0.58 





19a' (Casp2) 

- 

22a' 

0.26 




Table X. Calculated Spectrum of C-^oH-y-b’ 1 ' (^A") 


state 

E{kK) 

f 

polar 

3 3 A" 

16.2 

0.095 

X 

4 3 A" 

21.7 

0.030 

y 

1 3 A' 

24.2 

0.002 

z 

5 3 A" 

25.1 

0.049 

y 

6 3 A" 

26.2 

0.001 

xy 

7 3 A" 

29.2 

0.001 

xy 

8 3 A" 

30.2 

0.001 

xy 

10 3 A" 

34.6 

0.036 

X 

11 3 A" 

35.9 

0.047 

X 

12 3 A" 

37.4 

0.433 

X 

2 3 A 1 

28.4 

0.003 

z 

4 3 A 1 

34.0 

0.004 

z 


excitation 


3a“ (pi3) 

- 

5a" (pi5) 

5a“ (pi5) 

- 

6a" (pi6) 

3a" (pi3) 


8a" (pi8) 

19a' (Cbsp2) 

- 

6a" (pi6) 

4a" (pi4) 

- 

19a’ (Cbsp2) 

19a' (Cbsp2) 


7a" (pi7) 

19a' (Cbsp2) 

- 

9a" (pi9) 

19a' (Cbsp2) 

- 

8a" (pi8) 

2a" (pi2) 

- 

5a" (pi5) 

4a“ (pi4 ) 

- 

7a" (pi7) 

4a" (pi4) 

- 

7a" (pi7) 

4a" (pi4) 


6a" (pi6) 

3a" (pi3) 

” 

8a" (pi8) 

4a" (pi4) 

- 

6a" (pi6) 

3a" (pi3) 


6a" (pi6) 

5a" (pi5) 

- 

8a" (pi8) 

3a" (pi3) 

- 

6a" (pi6) 

5a" (pi5) 

- 

7a" (pi7) 

4a“ (pi4) 

- 

6a" (pi6) 

4a" (pi4) 

- 

8a“ (pi8) 

3a" (pi3 ) 

- 

7a" (pi7) 

2a" (pi2) 

- 

6a" (pi6) 

4a“ (pi4) 

- 

6a" (pi6) 

5a" (pi5) 

- 

7a" (pi7) 

4a" (pi4) 

- 

19a' (Cbsp2) 

19a' (Cbsp2) 

- 

6a" (pi6) 

2a" (pi2) 


19a 1 (Cbsp2 ) 

19a' (Cbsp2) 


7a“ (pi7) 

19a' (Cbsp2) 

- 

9a" (pi9) 

19a' (Cbsp2) 

- 

6a" (pi6) 


coeff . 

0.89 

0.72 

0.11 

0.40 

0.16 

0.11 

0.10 

0.10 

0.59 

0.34 

0.63 

0.11 

0.10 

0.54 

0.28 

0.57 

0.20 

0.54 

0.40 

0.49 

0.20 

0.11 

0.55 

0.19 

0.54 

0.18 

0.15 

0.51 

0.16 

0.14 


Table XI. Calculated Spectrum of C 10 H 7 -b* ( 2 A') 


state 

E(kK) 

f 

polar' n excitation 

coef f 

2 2 A" 

28.5 

0.006 

z 

19a' (Cbsp2) 

_ 

6a“ (pi6) 

0.38 





19a' (Cbsp2) 

- 

7a" (pi7) 

0.16 





5a" (pi5) 

- 

19a' (Cbsp2) 

0.13 





19a' (Cbsp2) 

- 

9a" (pi9) 

0.11 





19a’ (Cbsp2> 

- 

8a “ (pi8) 

0.11 

5 2 A ' 

29.5 

0.001 

y 

4a" (pi4) 

- 

7a" (pi7 ) 

0.66 





5a" (pi5) 

- 

6a" (pi6) 

0.14 

6 2 A' 

32.4 

0.005 

xy 

5a n (pi5) 

- 

7a" (pi7) 

0.41 





4a“ (pi4) 

- 

6a" (pi6) 

0.35 

l 2 A' 

33.7 

0.005 

xy 

4a" (pi4) 

- 

6a" (pi6) 

0.54 





5a” (pi5) 

- 

7a" (pi7) 

0.28 

8 2 A' 

34.5 

0.009 

y 

5a" (pi 5) 

- 

7a" (pi7) 

0.28 





5a" (pi5) 

- 

9a“ (pi9) 

0.22 





4a" (pi4) 

- 

8a* (pi8) 

0.18 





2a" (pi2) 

- 

6a" (pi6) 

0.11 

4 2 A" 

34.9 

0.002 

z 

19a’ (Cbsp2) 

- 

7a“ (pi7) 

0.44 





19a' (Cbsp2) 

- 

6a" (pi6) 

0.24 





19a' (Cbsp2) 

- 

9a" (pi9) 

0.11 

9 2 A ' 

37.0 

0.161 

y 

5a" (pi5) 

- 

6a" (pi6) 

0.84 


ib 



Table XII. Calculated Spectrum of C^QH 7 -b (■'■A') 


state 

E(kK) 

f 

polar 'n excitation 

coef f 

l 1 A" 

9.6 

0.001 

z 

19a' (Cbsp2) 

_ 

6a" <pi6) 

0.66 





19a' (Cbsp2) 

- 

8a” (pi8) 

0.26 

2^A" 

14.7 

0.004 

z 

19a 1 (Cbsp2) 

- 

7a" (pi7) 

0.60 





19a’ (Cbsp2) 

- 

9a" (pi9) 

0.20 





19a' (Cbsp2) 

- 

10a- (pilO) 

0.16 

3^' 

25.1 

0.003 

z 

19a* (Cbsp2) 

- 

8a- (pi8) 

0.62 





19a' (Cbsp2) 

- 

6a" (pi6) 

0.30 

2 X A ' 

30.0 

0.100 

xy 

5a" (pi5) 


6a" (pi6) 

0.78 

4 X A" 

30.6 

0.012 

z 

19a' (Cbsp2) 

- 

10a- (pilO) 

0.51 





19a' (Cbsp2) 

- 

7a“ (pi7) 

0.32 





19a' (Cbsp2) 

- 

8a" (pi8) 

0.10 

3^' 

33.1 

0.015 

xy 

5a" (pi5) 

- 

7a" (pi7) 

0.47 





4a" (pi4) 

- 

6a“ (pi6) 

0.33 





5a" (pi5) 


6a" (pi6) 

0.12 

4^' 

35.4 

0.045 

X 

19a' (Cbsp2) 

~ 

20a’ 

0.41 





19a' (Cbsp2) 

- 

21a' 

0.26 


*7 




Table XIII. Calculated Spectrum of C 10 H 9 _a+ < 1a ') 


state 

E(kK) 

f 

polar 'n excitation 

coef f 

2 1 A‘ 

25.8 

0.077 

xy 5a" (pi4 ) 

- 6a" (pi5) 

0.97 

3 1 A‘ 

29.0 

0.489 

x 4a"(pi3) 

- 6a" (pi5) 

0.95 


IS 


Table XIV. Calculated Spectrum of C 10 H g -a* ( 2 A") 


state 

E(kK) 

f 

polar 1 n excitat ion 

coef f 

2 2 A" 

15.7 

0.001 

X 

6a” (pi 5) 

— 

7a" (pi6) 

0.43 





4a" (pi3 ) 

- 

6a" (pi5) 

0.16 





5a" (pi4) 


8a" (pi7) 

0.10 

3 2 A" 

23.4 

0.006 

xy 

6a* (pi5) 

- 

8a" (pi7) 

0.29 





6a" (pi5) 

- 

9a" (pi8) 

0.19 





4a" (pi3) 

- 

8a" (pi7) 

0.10 





5a" (pi4) 

- 

6a" (pi5) 

0.10 

4 2 A" 

23.6 

0.002 

xy 

6a" (pi5) 

- 

9a" (pi8) 

0.26 





6a“ (pi5) 


8a" (pi7) 

0.19 





4a" (pi3) 

- 

7a" (pi6) 

0.10 

5 2 A" 

30.8 

0.014 

xy 

5a" (pi4) 

- 

8a" (pi7) 

0.49 





5a" (pi4) 


6a" (pi5) 

0.13 

1 2 A' 

34.0 

0.001 

z 

6a" (pi 5 ) 

- 

20a 1 

0.77 

6 2 A” 

34.6 

0.027 

y 

5a" (pi4) 

- 

6a" (pi5) 

0.25 





5a" (pi4) 

- 

7a" (pi6) 

0.17 





5a" (pi4) 

* 

8a" (pi7) 

0.16 





6a" (pi5) 

- 

8a" (pi7) 

0.14 

7 2 A" 

36.4 

0.239 

X 

4a" (pi3 ) 

- 

6a" (pi5) 

0.44 





6a" (pi5) 

- 

7a" (pi6) 

0.27 

2 2 A 1 

39.4 

0.002 

z 

6a" (pi5) 

- 

21a’ 

0.53 





6a" (pi5) 

- 

23a* 

0.21 





4a" (pi3) 

- 

20a ‘ 

0.10 



Table XV. Calculated Spectrum of C^gHg-a {^A') 


state 

E(kK) 

f 

polar 1 

'n excitation 

coef f . 

to 

> 

20.0 

0.025 

xy 

6a” (pi 5) 

- 7a" (pi6) 

0.57 





6a" (pi5) 

- 8a" (pi7) 

0.39 

3^' 

22.1 

0.417 

X 

6a“ (pi5) 

- 8a" (pi7) 

0.54 





6a" (pi5) 

- 7a H (pi6) 

0.40 

4^' 

34.9 

0.172 

y 

6a" (pi5) 

- 9a" (pi8) 

0.91 

3^" 

39.3 

0.002 

z 

6a" (pi5) 

- 22a 1 

0.47 





6a" (pi5) 

- 24a' 

0.31 





6a" (pi5) 

- 21a' 

0.18 


2 
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Table XVI. Calculated Spectrum of C 1 QH g -b + ( ^A' ) 


state 

E(kK) 

f 

polar 'n excitation 

coef f 

2 X A 1 

23.4 

0.184 

xy 

5a H (pi4) 

- 6a" (pi5) 

0.97 

3^' 

35.6 

0.138 

xy 

4a" (pi3) 
5a 11 (pi4) 

- 6a"(pi5) 

- 7a"(pi6) 

0.69 

0.26 

4 X A‘ 

40.0 

0.014 

xy 

3a" (pi2) 
5a” (pi4) 

- 6a" (pi5) 

- 8a" (pi7 ) 

0.73 

0.13 





Table XVII. Calculated Spectrum of C^QHg-b* ( 2 A“) 


state 

E(kK) 

f 

polar' n excitation 

coeff 

2 2 A“ 

13.3 

0.001 

xy 

5a“ (pi4) 

- 7a" (pi6) 

0.32 





6a" (pi5) 

- 7a" (pi6 ) 

0.26 

3 2 A" 

18.4 

0.006 

xy 

6a" (pi5) 

- 8a"(pi7) 

0.30 





4a" (pi3 ) 

- 8a" (pi7) 

0.29 





4a” (pi3) 

- 6a" (pi5) 

0.11 

4 2 A" 

22.8 

0.001 

y 

6a" (pi5) 

- 9a" (pi8) 

0.19 





5a“ (pi4) 

- 9a" (pi8) 

0.15 





6a" (pi5) 

- 7a" (pi6) 

0.12 

5 2 A" 

28.5 

0.087 

xy 

5a" (pi4) 

- 8a"(pi7) 

0.37 





5a“ (pi4) 

- 6a"(pi5) 

0.29 





6a" (pi 5) 

- 7a* (pi6) 

0.10 

6 2 A u 

30.6 

0.040 

xy 

5a" (pi4) 

- 7a" (pi6) 

0.33 





6a" (pi5) 

- 7a" (pi6) 

0.12 





4a" (pi3) 

- 7a“ (pi6) 

0.11 





5a" (pi4) 

- 9a" (pi8) 

0.11 

1 2 A' 

33.7 

0.004 

z 

6a* (pi5) 

- 20a' 

0.70 

7 2 A“ 

35.0 

0.005 

xy 

6a" (pi5) 

- 9a" (pi9) 

0.21 





4a" (pi3) 

- 9a" (pi9) 

0.15 





6a" (pi5) 

- 9a" (pi8 ) 

0.11 

8 2 A" 

38.8 

0.030 

y 

5a" (pi4 ) 

- 8a" (pi7 ) 

0.25 





5a“ (pi4) 

- 7a" (pi6) 

0.20 





4a" (pi3) 

- 6a* (pi5) 

0.15 





Table XVIII. Calculated Spectrum of C-j^Hg-b ( 1 A') 


state 

E(kK) 

f 

polar ' n excitation 

coef f 

2 X A‘ 

16.4 

0.108 

xy 

6a" (pi5) 

- 7a" (pi6) 

0.98 

3 L A' 

26.8 

0.037 

xy 

6a" (pi5) 

- 9a"{pi8) 

0.67 





6a" (pi5) 

- 8a"(pi7) 

0.30 

4 X A 1 

29.9 

0.375 

xy 

6a" (pi5) 

- 8a" (pi7) 

0.61 





6a" (pi5) 

- 9a"{pi8) 

0.31 

2 X A" 

37.0 

0.003 

2 

6a" (pi5) 

- 21a' 

0.84 
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Table XIX. Calculated Energies and Oscillator Strengths of Electronic 
Transitions of Naphthalene and Derivatives Relevant to the DIB Issue 


Energy (cm- 1 ) 

Wavelength 

(A) Oscillator 
Strength 

Species 

15,200 

6579 

0.105 

CioHe + 

20,500 

4878 

0.040 

+ 

CO 

X 

o 

o 

9,700 

10,309 

0.066 

C10H8' 

13,500 

7,407 

0.044 

CioH 8 ' 

18,600 

5376 

0.071 

CioH 8 ' 

16,100 

6211 

0.094 

CioH 7 -« + 

21,900 

4566 

0.030 

CioHya + 

9900 

10,101 

0.004 

C10H7-C*' 

16,200 

6173 

0.095 

C,oH 7 -P + 

21,700 

4608 

0.030 

C,oH 7 - 0 + 

9600 

10,417 

0.001 

C10H7-P- 

14,700 

6803 

0.004 

C10H7-P' 

15,700 

6369 

0.001 

CioHg-a 

20,000 

5000 

0.025 

CioHg-a' 

22,100 

4525 

0.417 

CioHg-a' 

13,300 

7519 

0.001 

CioHg-fl 

18,400 

5435 

0.006 

CioHg-0 

22,800 

4836 

0.001 

CioHg -|3 

16,400 

6098 

0.108 

CioHg-0" 





































